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ABSTRACT 

The known population of pulsars contains objects with four and five component pro¬ 
files, for which the peak-to-peak separations between the inner and outer components 
can be measured. These Q and M type profiles can be interpreted as a result of sightline 
cut through a nested cone beam, or through a set of azimuthal fan beams. We show 
that the ratio Rw of the components’ separations provides a useful measure of the 
beam shape, which is mostly independent of parameters that determine the beam scale 
and complicate interpretation of simpler profiles. In particular, the method does not 
depend on the emission altitude and the dipole tilt distribution. The different struc¬ 
tures of the radio beam imply manifestly different statistical distributions of Rwi with 
the conal model being several orders of magnitude less consistent with data than the 
fan beam model. To bring the conal model into consistency with data, strong effects of 
observational selection need to be called for, with 80% of Q and M profiles assumed to 
be undetected because of intrinsic blending effects. It is concluded that the statistical 
properties of Q and M profiles are more consistent with the fan-shaped beams, than 
with the traditional nested cone geometry. 

Key words: pulsars: general - pulsars: individual: J0631+1036 - Radiation mecha¬ 
nisms: non-thermal. 


1 INTRODUCTION 

In spite of a large and increasing number of detected ra¬ 
dio pulse profiles, a generic shape of pulsar beam remains a 
subject of debate. The mainstream models seem to support 
patchy or conal geometry. The patchy form is supported by 
the diversity and asymmetry of profiles, as well as by the in¬ 
voked distribution of individual components within the polar 
tube (Lyne & Manchester 1988; Manchester 2012). 

In a series of papers, (Rankin 1988; 1990; 1993, hereafter 
R93), Joanna Rankin provides arguments for approximate 
beam geometry in the form of two nested cones with an axial 
core component. It has been suggested that this beam geom¬ 
etry is approximately universal, with many pulsars having 
either the inner or outer cone, with a possible core compo¬ 
nent. A small group of profiles with 4 and 5 components (Q 
and M type, respectively) has been interpreted as a cut of 
sightline through both cones. The central component in M- 
type profiles is created by the additional co-axial core beam. 
In the works of J. Rankin, the angular radii of the cones have 
been estimated for pulse longitudes measured at the outer 
50% peak flux of components. At 1 GHz they are equal to: 
pfn = 4.3 °p~ 1 / 2 an d p®jj t = 5.8°P -1 / 2 for the inner and 
outer cone, respectively; (P denotes a pulsar period, which 
in these equations should be specified in seconds). The re¬ 


sult has been confirmed by other groups, who measured the 
conal pair widths at a different flux level and frequency (Gil 
et al. 1993, hereafter G93; Kramer et al. 1994, hereafter K94; 
Mitra & Deshpande 1999, hereafter MD). MD have tenta¬ 
tively identified three cones, out of which we select the inner 
two, because they outnumber the third-cone case, and their 
ratio is consistent with that derived in other studies. 

IWrightl (120031) introduced a special-relativistic model of 
drifting pulsar beams in which the two cones are associ¬ 
ated with two particular sets of dipolar magnetic field lines: 
the last open lines and the critical lines, which separate 
zones of opposite-sign charge at the light cylinder (of radius 
P lc = cP /(2tv), where c is the vacuum speed of light). In 
pulsar magnetospheres, these lines form two co-axial tube¬ 
shaped surfaces with different opening angles, 2pi 0 and 2p cr t. 
In a dipolar field with magnetic moment parallel to the ro¬ 
tation axis, it holds that Pcr(r)/pi 0 (r) = 0.75. If the ra¬ 
dial distance r from the center of a neutron star is not too 
large (r -C R\ c ), the result is independent of r. For a large 
tilt of stat ic-shaped di pole, the theoretical ratio increases 
up to 0.82. IWrightl f2003|) notes that the numbers are close 
to the size ratio of the cones observed by Rankin (1993): 
Rp° = pf® /pout ~ 0.74. In Tablc[l]we give other values of the 
ratio, as determined from observations by several research 
groups. A narrowing of cones with increasing frequency v, 
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MD 

R93 

G93 


K94 (P~ 

■“) 

K94 (P- 

- 0 . 5 ) 

V [GHz] 

1.0 

1.0 

1.4 

1.4 

4.75 

10.55 

1.4 

4.75 

10.55 

Rf [%} 

100 

50 

10 

10 

10 

10 

10 

10 

10 

Pin [°] 

4.1 

4.3 

4.9 

5.3 

4.5 

4.77 

4.9 

4.4 

4.5 

pout [°] 

5.1 

5.8 

6.3 

6.23 

5.76 

5.48 

6.3 

5.9 

5.5 

R p 

0.8 

0.74 

0.78 

0.85 

0.78 

0.87 

0.78 

0.75 

0.82 


Table 1. Angular radius of the inner (p; n ) and outer (p 0 ut) cone, and their ratio R p = Pi n /p 0 ut, as determined by various statistical 
studies of pulsar profiles. The values are based on profile widths measured at a different fraction Rp of the components’ peak flux. In 
the case of MD, who have identified three cones, we provide values for the inner two. The last three columns (from K94) and the values 
from MD refer to a fit with a fixed period dependence of 0,5 . 


can be inferred from the last 6 columns of Tab. [l] In spite 
of this, the ratio of cones’ size R p remains ^-independent. 
This is consistent with the cones occupying the same mag¬ 
netic held lines at different altitudes. The geometry of the 
emission region seems to follow the haring geometry of the 
dipolar magnetic held. 

The cone size ratio R p depends on the Hux level at which 
locations of components in a prohle are measured. This is 
usually set as a fraction Rf of peak hux of a considered com¬ 
ponent (see Fig. 4 in Kramer et al. 1994). For similar compo¬ 
nents of the inner and outer cone (similar width and shape), 
R p should slightly increase with decreasing Rf. Comparison 
of results from R93, G93 and K94 suggests a weak increase 
of R p for Rf decreasing from 50 to 10%. Columns 2-4, how¬ 
ever, which also include the result of MD for Rf = 100%, 
provide little evidence for this. In this paper we deal with 
the multicomponent prohles of class Q and M in which com¬ 
ponents often partially overlap with each other. We hnd that 
peaks of such blended components can, on average, be more 
easily identified than the points corresponding to a lower 
hux fraction. Therefore, to minimise the blending problems, 
we assume Rf = 100%. 

In addition to the patchy and conal beams, a variety 
of more complicated shapes have been considered, such as 
the hourglass shape (Weisberg & Taylor 2002), elliptic (MD; 
Perera et al. 2010), various systems of fan beams, eg. spoke¬ 
like or wedge-like (Dyks et al. 2010, hereafter DRD10; Wang 
et al. 2014; Teixeira et al. 2016) and in the form of a mod¬ 
elled polar cap rim as determined by the magnetic heldline 
tangency condition at the light cylinder (Biggs 1990; Dyks 
& Harding 2004). A hybrid form consisting of ‘patchy cones’ 
has also been considered by Karastergiou & Johnston (2007) 
and shown to reproduce some statistical properties of pulsar 
prohle ensamble. 

There is a growing evidence that pulsar radio beams 
generally do not have a conal geometry. DRD10 proposed a 
radio emission beam in the form of multiple fan beams cre¬ 
ated by plasma streams diverging from the magnetic dipole 
axis. This radio emission geometry, resembling the pattern 
of spokes in a wheel when viewed down the dipole axis, 
has shown many advantages when compared to the nested- 
cone case (see Dyks and Rudak 2012; Desvignes et al. 2013; 
Wang et al. 2014; Dyks & Rudak 2015). The new model has 
managed to explain the main frequency-dependent features 
of multicomponent prohles, such as the radius-to-frequency 
mapping and the relativistic core lag (Gangadhara & Gupta 
2001). It also provides a successful interpretation of double 


notches observed in some averaged prohles (McLaughlin & 
Rankin 2004). 

With the long-term monitoring, pulsar beams can be 
mapped for precessing objects, especially those which under¬ 
gone the fast precession caused by the relativistic spin-orbit 
coupling (Kramer 1998; Hotan et al. 2005; Clifton & Weis¬ 
berg 2008). The rare examples that have been mapped so far 
show that there is much to learn about the form of pulsar 
beams. Beam maps in Desvignes et al. (2013) and Manch¬ 
ester et al. (2010) suggest elongated patterns that point at 
the magnetic pole, in line with the fan-shaped pattern dis¬ 
cussed in DRD10 (hg. 18 therein). 

The majority of previous works have focused on a study 
of widths of radio prohles. The widths, however, are sensi¬ 
tive to a number of factors, such as the rotation period P, 
physical size of the emission region, altitude of emission, and 
the macroscopic pulsar geometry (dipole tilt and viewing an¬ 
gle). All of them are convolved and, except from the period, 
unknown for the majority of objects. 

The Q and M prohles provide a useful tool for decipher¬ 
ing the pulsar beam shape, because it is possible to study 
the ratio of separations between their components, instead 
of the full width of a prohle. With a good accuracy, distri¬ 
butions of such ratios are insensitive to major parameters 
that complicate the analysis of prohle widths (such as the 
emission altitude, frequency v, rotation period P, and the 
specihc form of a dipole tilt distribution). In this paper we 
study the statistics of components’ separation ratio for the 
conal and fan-beam models, and compare with the observed 
distribution. 

The outline of this paper is the following: In Section [2] 
we describe the main idea of the paper, and apply it to the 
conal model. Section [3] presents the way in which the ob¬ 
served distribution of the peak-separation ratio was deter¬ 
mined. In Section[4]we compare the theoretical distribution 
of the peak-separation ratio in the nested cone model with 
the observations. Section [S] does the same for the stream-like 
geometry, and is followed by discussion and conclusions. 


2 PEAK-SEPARATION RATIO IN THE 
NESTED CONE MODEL 

The idea is based on the measurement of the peak-to- 
peak width Win and Wont, for the inner and outer pair 
of conal components, respectively, to derive the width ra¬ 
tio Rw = Win/Wont- For a central cut through the nested 
cones with angular radii pin and p ou t, the ratio Rw is max- 
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Figure 1. Top half of a nested cone beam (two half-circles) with 
the cone size ratio R p = 0.75. The set of horizontal lines (paths of 
different sightlines) presents intervals of viewing angle that corre¬ 
spond to a fixed interval of the peak-separation ratio A Rw =0.1. 
Values of the impact angle /3 (in units of the outer conal radius 
Pout) are given on the left, values of Rw - on the right hand side. 
Note that it is a lot more probable to observe R\y S (0.65, 0.75) 
than the smaller values. 


imal and equal to the ratio of cones’ size: R p = pin/pout• 
This is also the most likely value, which should be vastly 
dominant in the data. Fig. |T] presents the upper half of a 
nested cone beam with R p = 0.75, viewed down the dipole 
axis. The horizontal lines mark the paths of the line of sight 
for observers located at different impact angles (3 = C, — a, 
where C, is the angle between the sightline and rotation axis 
fi, and a is the tilt of a magnetic dipole with respect to 
fi. The real spherical geometry of the sightline cut implies 
that the viewing paths are curved in general. However, the 
straight lines of Fig. [T| provide a good approximation when¬ 
ever a pout and p ou t <C 1 rad. The paths have been 
plotted for the equidistant values of the peak separation ra¬ 
tio Rw, printed on the right-hand side. Adjacent values of 
of /3/pout (shown on the left) quickly approach each other 
with decreasing Rw , which means that the chance to observe 
the inner components at a small separation is considerably 
smaller than observing them at a larger distance. For exam¬ 
ple, the probability to observe Rw in the range (0.25, 0.35) 
is (0.73 — 0.71)/0.48 = 0.04 times smaller than to observe a 
larger Rw in the same-width interval of (0.65,0.75). 

For a nested cone beam with a universal R p , the dis¬ 
tribution of Rw is a single-peaked function, monotonously 
increasing towards the sharp peak at Rw = R P - As demon¬ 
strated in Appendix m the function n(Rw) can be easily 
derived in the flat case shown in Fig. [l] 


n(Rw) = C 


dn 

dRw 


Rwp-Rf) 

Upo »\Rt-R* w y/2(l-Rl w )Z/2 


(1) 


where C is a normalisation constant. 

This distribution has important advantages over a di¬ 
rect statistical study of pulse widths, because Win and Wout 
are expected to depend on emission altitude, observation 
frequency, rotation period of the star, and the dipole tilt a. 
In the case of small beams, typical of ordinary pulsars, the 
flat case of eq. 0 presents a good approximation within 
most of the (a, £) parameter space, except from cases of 
the nearly aligned geometry. The latter, however, are not 
numerous, so the actual distribution n(Rw), as calculated 
with the strict account of the spherical trigonometry, is close 
to eq. 0. Therefore, it is worth to discuss the properties of 



R w 


Figure 2. Dependence of the Rw distribution on period, with a 
and £ sampled isotropically. Note that the grey histogram (P = 1 
s) essentially coincides with the analytical case of eq. o (p = oo, 
marked with the plus signs), despite the former includes cases 
with ol < 90°. The solid line case for P = 0.1 s is very close to 
the aforementioned (long P) cases. All of these distributions are 
therefore appropriate for majority of normal pulsars. The num¬ 
bers on the vertical axis refer to P = 1 s, with the other his¬ 
tograms normalised at the same peak value. 


the Rw distribution within the range of validity of eq. ©. 
If the cones are associated with the last open and critical 
field lines (or any lines defined by two fixed values of the 
footprint parameteiQ), then the value of Rw is not expected 
to depend on the rotation period P or the emission altitude 
r. This is because a change in P or r just rescales the beam 
shown in Fig. |T] with no influence on the relative geometry 
of cones and the statistics of Rw- For the aforementioned 
B-field lines, the ratio pi n /pout does not depend on r as long 
as r Rio- Therefore, we fix it at R p = 0.75. A choice of 
frequency v should not affect Rw either, if the peak emis¬ 
sion at different v corresponds to different altitudes, but the 
same type of dipolar field lines (last open/critical). Note 
that Pout(P,r) in eq.0 can be incorporated into the nor¬ 
malisation constant C. This is because the fraction in eq. 0 
(let us denote it by F, so that n = CpoutF) is almost in¬ 
sensitive to P and r. A ‘total’ distribution that incorpo¬ 
rates beams of different size p out ,i can then be written as 
n « J2i Cpout.iP = CFy\ Pout,i = CF. For all the afore- 
described reasons, the distribution of peak separation ratio 
Rw provides a useful, one-dimensional tool for testing the 
pulsar beam shapes. It is insensitive (or very weakly sensi¬ 
tive) to the uncertain parameters, and allows us to avoid the 
usual two-dimensional analysis (eg. of n{W,P)). On the bad 
side, it is applicable only for a limited number of objects 
with four and more components (of Q and M type). 

Since Fig.Q]and eq. 0 are only valid for large a S> p ou t, 
we have determined the distribution n(Rw) numerically, by 
calculating exact values of Rw for a large sample of beams 


1 The footprint parameter is a ratio of the magnetic colatitude of 
an arbitrary point and the magnetic colatitude of the open field 
line boundary, measured at that point’s radial distance. 
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(10 7 ) simulated for isotropically distributed angles of a, and 
The opening angle of the outer cone has been calculated 
with the usual dipolar formula: 

cospout = (2 - 3r/f?i c )(4 - 3r/7?ic) _1/2 . (2) 

A variety of periods P, typical of normal pulsars has been 
tried to verify the near-independence of n(Rw) on P (see 
Fig.®. Strong discrepancy has appeared only for P less than 
a few tens of milliseconds (the dotted line in Fig. [5] is for 
P = 10 ms), whereas the periods in our sample of observed 
Q and M pulsars range between 0.16 and 2.25 s, with an 
average of 0.84 s. For smaller P, the pulsar beam is larger, 
and the curvature of the sightline’s path within the beam 
is more important. As shown in Appendix m this effect is 
second order in p ou t, ie. it is usually small. Because of the 
curved viewing paths, Rw decreases, hence the Rw values 
from the highest histogram bin (Rw £ (0.7, 0.75)) start to 
pour over to adjacent bins on the left (with Rw < 0.7). See 
the dotted line in Fig. ® This effect is more pronouced for 
narrower histogram bins. 

Following R93, the value of r has been set to 220 km, 
and we have assumed pi n = 0.75p ou t- The period has been 
set to P = 1 s, which is a round number close to the average 
P in the observed sample of Q and M pulsars. The average P 
of the total population of known pulsars is smaller, however, 
evidence has been presented for that the Q and M profiles are 
mostly observed in old objects (Rankin 1990). The resulting 
pulse width for the inner and outer cone has been calculated 
with the spherical cosine theorem 

cos(W;/2) = (cos pi — cosacosC)(sinasin<J) _1 , (3) 


where the index i refers to ‘in’ or ‘out’. Since a and £ are 
blindly sampled from an isotropic distribution, for most of 
them the line of sight does not traverse both cones. This 
happens only when eq. m gives finite solutions for both 
Wi, thus we accept only those (a,£)-pairs for which 

\P\ < P'\n and cos p ou t > cos (a + £). (4) 


For P = Is and r = 220 km, pi n = 4.38°, so the conditions 
are passed by just a few percent of total number of cases. 

Fig. 0 presents the n(Rw) distribution for selected val¬ 
ues of dipole inclination a and for two values of P = 1 s (top) 
and 0.1 s (bottom). In the long-period case (top panel), even 
moderate dipole inclinations result in a distribution which is 
well described by eq. ©: the solid line histogram for a = 45° 
is indiscernible from the grey orthogonal (a = 90°), or from 
the analytical case). Pronounced difference can only be seen 
for a very small inclination a = 10° and short P (Fig. ®:>, 
dotted). The reason for this is that for decreasing a both Win 
and Wont scale approximately as Wi / sin a which makes Rw 
very stable. In Appendix[B]a second-order expansion of Rw 
is made for the special case of the central cut (a = £): 


Rw 


R P 


1 - K pint \ 
24 sin 2 a J 


= Roll- 0.018 


2 

Pout 


,(5) 


where the numeric value on the right corresponds to R p = 
0.75. As one can see, Rw starts to perceptibly depend on a 
only for a nearly-aligned geometry (a < p 0 u t). According to 
eq. El when a is decreasing, Rw decreases with respect to 
the flat-case value of R p = 0.75. Hence Win is increasing a 
bit slower than W ou t- Because of the square dependence on 
pi/ sin a, a considerable divergence may appear only when 



Figure 3. Dependence of the Rw distribution with the dipole 
tilt a for P = Is (top) and P = 0.1 s (bottom). Only the viewing 
angle £ was sampled isotropically. Note that the Rw distribution 
practically does not depend on a, except from when both a = 10° 
and P = 0.1 are simultaneously small (dotted line in bottom 
panel) which is a rare circumstance. 


pi is relatively large or sin a small. In the rare cases when 
sin a p 0 nt (dotted line in Fig. Eb), a strong discrepancy 
from the flat analytical case appears and the histogram does 
not extend all the way up to R p . In the case of such a small 
a the line of sight is capable of staying for most of the time 
between the cones. An extreme example with Win ~ 0 and 
Wout « 27 t, is the case with the H axis located half way 
between the cones (a = pin + (p 0 nt — Pin) /2) and ( « (pout — 
pin)/ 2. However, since such cases are rare they do not affect 
the total Rw distribution which is based on the isotropic 
distribution of a. Therefore, the sensitivity of n(Rw) to the 
key geometrical parameter, the dipole tilt a, is marginal, and 
nearly completely reduced as compared to the sensitivity of 
the n(W) distribution. 


3 THE OBSERVED REDISTRIBUTION 

To verify the simulated distribution of Rw we have com¬ 
puted the observed distribution (Fig.®:), by estimating the 
peak separations for 30 pulsars of Q and M class. The proce¬ 
dure started with a selection of as many Q and M-type pro¬ 
files as possible, followed by a visual estimate of component 
number and location. The objects have been selected from 
the works of R93, G93, Hankins and Rankin (2010), and 


2 Eq. l® has been derived for p Q ut 1 rad and sin a p Q ut • 


© 0000 RAS, MNRAS 000, 000-000 





























Peak separation ratio as a probe of pulsar beam 5 


Name 

Ref. 

Type 

Qual. 

P 

s 

Freq. 

GHz 

01, 02, 071 — 1 , 071 

o 

Win 

O 

Wout 

O 

Rw 

B0329+54 

1 

M 

i 

0.714 

1.41 

33.0, 41.3, 49.8, 54.2 

8.5 

21 

0.40 

J0401-7608 

2 

Q 

2 

0.545 

0.6 

-4.2, -0.7, 3.5, 7.5 

4.2 

11.7 

0.36 

B0621-04 

3 

Q 

1 

1.039 

1.408 

162.5, 168.0, 174.6, 179.8 

6.6 

17.3 

0.38 

J0631+1036 

4, 5, 6 

Q 

1 

0.287 

1.4 

-9.2, -2.8, 3.0, 8.5 

5.8 

17.7 

0.33 

J0742-2822 

7, 8 

M 

2 

0.166 

1.5 

-7.6, -6.1. -1.5, 2.5 

4.6 

10.1 

0.46 

J1034-3224 

4, 9 

Q 

2 

1.150 

0.6 

-33.0, 0.0, 14.0, 39.0 

14.0 

72 

0.19 

B1055-52 

4 

Q 

1 

0.197 

3 

-89.0, -79.0, -66.3, -59.0 

12.7 

30 

0.42 

B1237+25 

10 

M 

1 

1.382 

1.4 

42.45, 44.65, 50.7, 52.0 

6.0 

9.5 

0.63 

J1326-5859 

4 

Q 

2 

0.477 

3 

-7.9, -2.0, -0.2, 5.0 

1.8 

12.9 

0.14 

J1327-6222 

8 

Q 

2 

0.529 

3 

-2.5, 0.0, 3.2, 4.6 

3.2 

7.1 

0.45 

J1536-3602 

4 

Q 

2 

1.319 

0.6 

-6.5, 1.0, 7.9, 14.3 

6.9 

21 

0.33 

J1651-7642 

4 

Q 

2 

1.755 

1.5 

-1.0, 2.5, 6.7, 11.0 

4.2 

12.0 

0.35 

B1737+13 

3 

M 

2 

0.803 

1.408 

180.0, 183.5, 192.0, 195.0 

8.5 

15 

0.57 

B1738-08 

3 

Q 

1 

2.043 

0.61 

173.8, 177.0, 182.7, 185.8 

5.7 

12 

0.47 

B1804-08 

3 

Q 

1 

0.163 

1.408 

151.0, 158.7, 166.2, 169.3 

7.5 

18.3 

0.41 

J1819+1305 

11 

Q 

1 

1.060 

0.327 

-11.5, -5.3, 7.2, 11.5 

12.5 

23 

0.54 

B1821+05 

12 

M 

2 

0.752 

0.43/1.401 

-13.0, -6.0, 6.5, 10.5 

12.5 

23.5 

0.53 

B1831-04 

3 

M 

1 

0.290 

0.606 

138.5, 153.0, 212.5, 233.5 

59 

95 

0.63 

B1845-01 

3, 10 

Q 

2 

0.659 

1.42 

62.0, 64.2, 68.4, 74.4 

4.2 

12.4 

0.34 

B1857-26 

3 

M 

2 

0.612 

0.61 

169.0, 176.5, 192.3, 198.5 

15.8 

29 

0.53 

B1910+20 

12 

M 

2 

2.232 

1.615 

-7.0, -3.3, 4.5, 7.4 

7.8 

14.4 

0.54 

B1918+19 

13 

Q 

2 

0.821 

0.327 

-24.7, -10.7, 1.0, 23.5 

11.7 

48 

0.24 

J1921+2153 

4 

Q 

2 

1.337 

3 

-3.7, -1.8, 1.2, 3.1 

3.0 

6.8 

0.44 

B1929+10 

1. 14 

Q 

2 

0.226 

1.71 

205.5, 207.5, 210.5, 213.4 

3.0 

7.9 

0.38 

B1952+29 

3 

Q 

2 

0.426 

1.408 

253.6, 259.0, 264.7, 268.8 

5.7 

15.2 

0.37 

B2003-08 

3 

M 

2 

0.580 

0.408 

162.0, 174.0, 205.0, 214.0 

31 

52 

0.60 

B2028+22 

12 

Q 

1 

0.630 

0.430 

-4.8, -0.7, 2.6, 6.9 

3.3 

11.7 

0.28 

B2210+29 

3, 10 

M 

2 

1.004 

1.42 

160.3, 162.6, 171.8, 174.9 

9.2 

14.6 

0.63 

B2310+42 

1 

M 

2 

0.349 

1.41 

106.7, 109.0, 115.3, 117.8 

6.3 

11.1 

0.57 

B2319+60 

3 

Q 

2 

2.256 

0.925 

173.9, 179.0, 182.3, 188.8 

3.3 

14.9 

0.22 


Table 2. A subset of Q and M-type pulsars for which the peak separation ratio Ft\y was possible to estimate. The component phases 0j 
refer to an arbitrary zero phase. References: 1) von Hoensbroech &: Xilouris (1997); 2) Johnston et al. (1998); 3) Gould & Lyne (1998); 4) 
the ATNF database (http://www.atnf.csiro.au/people/joh414/ppdata); 5) Weltevrede et al. (2010); 6) Teixeira et al. (2016); 7) Johnston 
et al. (2005); 8) Karastergiou & Johnston (2006); 9) Manchester et al. (1998); 10) Seiradakis et al. (1995); 11) Rankin & Wright (2008); 
12) Hankins &; Rankin (2010); 13) Rankin et al. (2013). 


the other sources itemised in the caption to Tabled Their 
profiles were grouped into three classes of quality (mostly de¬ 
termined by the easiness to discern components), with the 
worst group rejected. Profiles of the remaining 30 objects, 
listed in Table [2] were viewed at a large (clear) scale, with 
four vertical lines overplotted at the suspected locations of 
components’ peaks. In five-component profiles the central 
one has been ignored. The coordinates of the vertical lines 
were used to calculate Win, W ou t and Rw = Wm/Wout given 
in Table O In some cases (like B1821+05) it was helpful to 
refer to two frequencies to resolve doubts about the existence 
or location of a specific component. 

The profiles have purposedly been not decomposed by 
fitting analytical functions, mainly because such functions 
are unknown, and different components are apparently de¬ 
scribed by different analytical shapes (cf. fig. 1 in Kramer et 
al. 1994 with the outer components of J0631+1036 in fig. 2 
of Weltevrede et al. 2010). Intrinsic intramagnetospheric ef¬ 
fects often appear to make the components asymmetric, and 
the drifting phenomenon can possibly make them roughly 
triangular or trapezoidal. A fitting of the Gauss curves can 
therefore give false results, biased by the use of the incor¬ 
rect function in the decomposition. In the case of blended 
asymmetric components, the unguided Gaussian fitting is 
incapable to provide precise estimate of the components’ 


positions or even of their number (the latter must be de¬ 
cided by eye also for isolated asymmetric components, see 
the fit of the von Misses functions in fig. 2 of Weltevrede 
& Johnston 2008). For these reasons we have identified the 
components visually. The rather large width of bins used 
in our observed Rw-histogram (ARw = 0.05) makes the 
analysis less sensitive to errors. Moreover, the difference in 
the predictions of the models will be shown to be so large, 
that even the approximate estimate of the observed Rw his¬ 
togram is useful. 


Our identification of the morphological type (Q or M) is 
based solely on the number of easily identifiable components, 
so it may be different from that of R93, who also considered 
spectral and circular polarisation properties. Out of the 17 
pulsars common for our Tab. [2] and R93, only 8 have been 
designated a definite type in R93. Three of those (B0621—04, 
B1845—01 and B1918+19) have a different type assigned 
(respectively M, cT, and cT in R93). If all the three objects 
are rejected from the analysis, the observed Rw distribution 
retains its boxy shape with little impact on our conclusions. 
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Figure 4. Comparison of the simulated Rw distributions (conal 
model in panel a; fan beam model in b) with the one observed 
for Q and M type pulsars (panel c). The numbers in a and b give 
the KS probablity of consistency with the observed distribution. 
Note the low likelihood of the conal model. Grey parts of the 
histograms in b and c denote the M-type profiles. An unknown 
fraction of these is also included in the histogram of panel a, which 
shows the sum of the Q and M profiles. 

4 CONAL MODEL VERSUS DATA 

The observed distribution is presented in Fig. whereas 
the one simulated for the nested-cone model - in Fig. [4]t. 
The distributions are completely different. The conal model 
distribution is dominated by the value Rw = R P , which 
corresponds to the beam size ratio. The observed n(Rw) 
peaks at the value of Rw ~ 0.375 which should have been 
nearly absent in the data. The Kolmogorov-Smirnov (KS) 
test (Press et al. 1992) excludes the common origin of the 
distributions, giving it a probablity of 1CP 12 . 

The total distribution of the conal model is marginally 
consistent (common origin prob. of 0.04) with the M-type 
part of the distribution alone (grey part of the histogram 
in Fig. 2};). However, we find no convincing reason to argue 
that the nested cone structure is only responsible for the ob¬ 
served M-type profiles, whereas the Q profiles have different 
origin. The ratio of Q to M pulsar numbers is difficult to es- 
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Figure 5. The influence of selection effects caused by component 
blending on the conal Rw distribution. The grey histograms are 
affected by the limited resolving capabilities: a) A0 r es = 0.2W ou t; 
b) A(f> Tes = 2.1°. The width of the affected histograms (or their 
humps) can be estimated from eqs. ([6]) and © (dashed and dot- 
dashed vertical line, respectively). The unaffected histogram from 
Fig. UK is shown as a reference. 

timate in the conal model, because the detectability of the 
core depends on the central beam parameters and telescope 
sensitivity. Therefore we focus on the total distributions. 

Since the theoretical distribution is dominated by the 
peak at R p , the observed n(Rw) could only be explained 
by the nested-cone beams if most of them has the size ratio 
R p ~ 0.3 — 0.65 with most common values within (0.3, 0.45). 
This is not consistent with the findings described in the in¬ 
troduction (Tab. [I]). Even in the case of MD, who tentatively 
identified additional large cone with p = 1.3p 0 ut, their data 
are dominated by the ratio R p = 0.8, inconsistent with the 
observed statistics of Rw ■ Furthermore, we are not aware of 
any magnetospheric arguments, which could support such a 
large variety of R p , as implied by Fig. 0K. 

4.1 Selection effects in the conal model 

Individual observed components may be difficult to distin¬ 
guish because of their considerable width, or a large noise 
(instrumental or intrinsic, ie. related to a large dynamic 
range of dissimilar single pulses). It is therefore reasonable 
to introduce a threshold for their minimum resolvable prox¬ 
imity A 4> le s- This can be defined either as a fraction ew of 
the total pulse width W ou t or as an absolute value measured 
in degrees. In the first case the Rw histogram becomes nar¬ 
rower and less spiky, as shown in Fig. [SJu The part of the 
histogram left of the vertical dashed line becomes empty be- 
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cause Win < ewWout (ie. Rw < ew) for the peripheral cut 
through the beam. Provided ew > (1 — R P )/ 2 = 0.125, the 
part on the right-hand side of the dot-dashed line also van¬ 
ishes because the distance between adjacent inner and outer 
components (fWout — Win)/2) is smaller than the resolution 
ewWout for the central cuts through the beam. The result of 
Fig. [S]i does not depend on a since the resolution is scaled 
just as the profile width. For ew = 1/3 the histogram con¬ 
tracts to the zero width (the dashed and dot-dashed lines 
meet at Rw = 1/3) and no components can be resolved in 
a profile. 

Assuming that we are incapable to resolve components 
at a distance smaller than ewW ou t, the conal distribution 
can be made consistent with observations for ew ~ 0.22 
(common origin KS probablity: 0.58; the probability stays 
above 10 -2 within 0.18 < ew < 0.25). The observed dis¬ 
tances between adjacent components (A(/ a< jj) indeed have 
a distribution (not shown) which decreases steeply below 
A^adj/Wout ~ 0.2. It may therefore be possible that the 
observed distribution of Rw (Fig. [3J;) is heavily distorted 
by our limited capability to resolve intrinsically overlapping 
components. Such profiles, however, (with the outer com¬ 
ponents merged with their inner neighbours), should either 
have two well-separated components (D class), or three com¬ 
ponents with the central one well-separated (by 0.375VF O ut) 
from the inner conals (class T). These should be numerous, 
because they represent the right hand-side peak in the conal 
histogram. This would imply that conal components in sev¬ 
eral profiles of D and T class should consist of two (inner 
and outer) blended components. 

In the case of the fixed (Wo U t-independent) resolving 
capability the shape of the Rw histogram (obtained for all 
combinations of isotropically distributed a and </) depends 
on the value of A</ rea (equal to 2.1° in Fig. [5}o). The reason 
can be readily seen by decomposing the total Rw histogram 
into subhistograms which correspond to a fixed value of a 
(or a narrow interval of a) and the isotropic For a small 
a, profiles are wider because they are viewed at small angles 
with respect to the rotation axis, hence, the resolving lim¬ 
itations apply only for narrow profiles observed at a larger 
£ and a. The resulting Rw histogram is then composed of 
the unaffected part (with the sharp peak at Rw = R P , cor¬ 
responding to the circumpolar viewing at a small </) and 
a range of narrower sub-histograms corresponding to larger 
viewing angles (affected by A(/ res ). The strongest (most nu¬ 
merous), and narrowest contribution comes from the equa¬ 
torial viewing (cases with £ ~ a ~ 90°) and is visible in 
Fig. [5)o as the protruding part between the dashed and dot- 
dashed lines. 

Locations of the resulting bumps in the histogram can 
easily be determined analytically. The dashed line constrains 
the region where Win + A</ res , which corresponds to 


oi 2 A </: 

P ^ Pin 


2 

res 


Rw 


A</i 
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(6) 


The dot-dashed line sets the upper limit for the region 
wherein (Wont - W / ' in )/2 ^ A </> lea , ie.: 
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At/res + Pin ~ P~o ut 
A (/res + Pout Pin ' 


(7) 

( 8 ) 


In the case of large A</ res and a ~ 90°, the limiting condi¬ 
tions © and © may exclude the entire parameter space (all 
profiles unresolved/rejected, A(/ a dj < A(/ res for any 0). This 
happens when A(/ 2 es = (p 2 ut — p 2 n )/2 and Rw = 1/3. The 
histogram then has a single bump (or break) at Rw = 1/3, 
and the two vertical lines in Fig. [5] coincide. The actual 
look of a distribution affected by the limited resolution de¬ 
pends on the relative value of A</ rea as compared to the 
scale of the beam. The width of the humps (or of the his¬ 
togram itself) will change whenever parameters such as the 
independent r, the rotation period, or the lateral bound¬ 
aries of the cones are changed. Tests performed for different 
values of the fixed A</ rea , have shown that the undistorted 
part of the histogram (with the peak at Rw = 0.75) usu¬ 
ally contributes considerably to the overall shape, and the 
probability of consistency can hardly exceed 6 x 10 -4 (at 
A(/ re s = 2.5° and P = 1 s). 

The conal model is then found a poor representation of 
data, unless a properly-tuned selection effect is called for: it 
should be impossible to resolve components separated by less 
than 0.2IT ou t, regardless of the profile width Wout- This is 
possible if the intrinsic width of components, responsible for 
the blending, corresponds to a fixed fraction of polar tube. 
However, for the best-fit value of ew = 0.22, the number of 
profiles with unresolved peripheric components (that would 
be classified as the apparent double D or triple T profiles) is 
four times larger than the total number of known (resolved) 
Q and M profiles. That would imply that for some ~ 200 
pulsars of D and T type, the outer components are composed 
of the unresolved pairs. 


5 PEAK SEPARATION RATIO IN THE 
STREAM MODEL 

To learn the shape of the Rw distribution for the stream 
model of pulsar beam (DRD10; Dyks & Rudak 2012; Wang 
et al. 2014) we consider the simple geometry of emission lim¬ 
ited to separate magnetic azimuths <j> m . Nearly all observed 
radio profiles have at most five components. Therefore, for 
each beam in the sample, ten radio-emitting streams, in a 
‘5+5’ fashion, is assigned to the polar region: five of them 
in the upper (poleward) half of the polar tube, and the 
other five in the lower (equatorward) part (see Fig. [©). The 
magnetic azimuths (f>m,i were selected randomly. For each 
azimuth, a uniform radio emission was assumed within a 
limited range of angles 6 m measured from the dipole axis: 
9™ in < dm < 0/ILx- As before, the values of a and / have 
been sampled isotropically. The analysis was limited to Q 
and M profiles, ie. we discarded all the cases with less than 
four intersections between the azimuth of the emitted beam 
and the sightline path. Cases with more than five cross¬ 
ings, which rarely appear for the extremely aligned geome¬ 
try (a ~ £ & #i7Lx), have also been ignored. The component 
separations Win and Wout have been calculated for the in¬ 
ner and outer pair of the crossing points, respectively (in the 
cases with five intersections, the central one, corresponding 
to the ‘core’ component, was ignored). The widths Wi have 
been calculated using the strict spherical trigonometric for¬ 
malism, as described in DRD (see eq. 19 therein). Some 
technical complications are also discussed in Appendix D of 
the present paper. 
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Figure 6. A view of a simple fan beam down the dipole axis fl. 
The thick radial sections emit radio waves within the range of 
magnetic angles 0 m , S (©(' ; ' n ' ^max)* The observer’s line of sight 
revolves around the star rotation axis fl, along the solid line path 
(arc with the arrow). The distance between S2 and /T corresponds 
to the dipole tilt a. Poleward viewing is presented (sin < sin a) 
of a pulsar with small a. In the case shown, a two-component 
profile of D class would be recorded. Note the definition of the 
magnetic azimuth (j> m . The dotted projections of magnetic field 
lines are separated by 10° in (i> rn ,. 

The Rw distribution calculated for the stream model is 
shown in Fig. [Tjj. It is different from the observed one (KS 
probability of consistency: 0.002), albeit it is nine orders of 
magnitude more probable than the raw conal distribution 
(by ‘raw’ we mean the distributions unaffected by the limited 
resolving capability). As in the observed case, the fraction 
of the M type objects (grey part) increases with Rw- This 
is because the extra space needed for the central component 
makes the leading side components (outer and inner) more 
distant from the trailing pair. Accordingly, Rw is closer to 
unity. The increased fraction of M-type profiles at large Rw 
is also expected for the conal model, because detection of the 
core requires a more central traverse through the beam. It 
needs to be emphasized, however, that the number of radio- 
emitting streams that exist in magnetospheres of different 
pulsars likely varies between 0 and ~5 (per one, poleward 
or equatorward, magnetic hemisphere). The cases with the 
small numbers of streams (1, 2, 3) are likely responsible for 
majority of the single, double and triple profiles. Therefore, 
the observed ratio of M and Q profiles, is also affected by 
the ratio of objects that actually have 4 or 5 streams, and 
not only by the statistics of the traverse through the beam 
with 5-streams0 

3 The multiparameter modelling of the relative numbers of dif- 



Figure 7. The influence of component blending on the Rw dis¬ 
tribution in the fan beam model. In a) the minimum distance of 
resolvable components is Abres = 0.2Wout- In b) Arj> Tes = 1.5°. 
Note that the limited resolution of blended components strongly 
increases the fraction of Q type profiles (white part of the his¬ 
tograms) . 


5.1 Selection effects in the stream model 

As before, two types of selection effects have been applied 
to the stream model: the incapability to resolve components 
located closer than ewWout and the IT-independent reso¬ 
lution limit of A (p rea (see Fig. Qr and b, respectively). The 
effect of these on the Rw distribution was similar to the 
one described for the conal model: in the case of the res¬ 
olution proportional to the width (fixed ew, Fig. [7^) the 
distribution becomes narrower, and subdistributions of Rw 
calculated for different a are the same (a-independent). At 
ew =0.11 the probability of consistency with data reaches 
0.04. 

For a fixed A(f> les the selection effects do not operate 
at small £ so the resulting distribution consists of several 
contributions of different widths (Fig. [TJo). The fixed- Arf> Tes 
distribution is narrower, however, probably on the account 
of its leftward skewness, it is not very consistent with the ob¬ 
served one (Fig. U) at any value of A0 reB (maximum prob¬ 
ability of consistency in the KS test: 0.016 at A</> res = 1.1°). 


ferent profiles (of S, D, T, Q, and M class), with the allowance for 
different numbers of streams in different objects, is a complicated 
subject which will be discussed elsewhere (Frankowski et al. 2016, 
in preparation; see also Karastergiou & Johnston 2007). 
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5.2 Influence of emission region geometry 

Contrary to the conal model (see Sect. [TJ, the geometry 
of the stream-shaped region of radio emission is not even 
weakly constrained, either by theory or observations. There¬ 
fore, we have probed parts of the available parameter space 
by varying the following parameters: 1) the minimum an¬ 
gular distance of the emitting streams from the dipole axis 
which has so far been set to zero); 2) a minimum az¬ 
imuthal distance of the streams in the magnetic azimuth 
d(j>m', 3) an interval A0 m of magnetic azimuths that are 
available for positioning the streams, ie. the original choice 
of the equatorward interval —90° < c j> m < 90°, and the 
poleward one 180° — 90° < <f) m < 180° + 90°, has been 
replaced with two narrower intervals which do not extend 
that far from the main meridian: —A <f> m < <j> m < A <f>m, and 
180° — A 4> m < 4>m < 180° + A0 m , with 0 < A <p m < 90°. 

The value of 9™ in was varied in the full range between 
0 and 9™^. With the increase of 9™ in the original distri¬ 
bution (shown in Fig. [4j>) transforms into one which peaks 
near Rw = 0 and decreases monotonically at larger Rw- 
The consistency with observations stays at the level of a 
fewxlO -3 , until 9™ in reaches 75% of 9™^. For a larger 9™ in 
the initially spoke-like shaped pattern of elongated streams 
starts to resemble ‘patchy cones’, and the probability of con¬ 
sistency with data drops down to 2 x 10~ 4 at 9™ in = 0.9 
For all results included in this paper, the value of 9™ ax was 
set equal to the conal value of p out , although we have also 
tried the low multiplicities ip 0 u t, with i between 1 and 5. 
Since 9™^ just rescales the beam, the result does not de¬ 
pend on 9m ax as long as 9^ x <C 1 rad and the selection 
effects are neglected. 

The limit for the minimum allowable azimuthal separa¬ 
tion of streams makes the agreement with data worse: the 
distribution tends towards a narrow bump at Rw =0.2. 

By keeping the streams closer to the main meridian 
(A<j>m = 0.75 x 90°), the probability of consistency with 
the observations can be increased to a considerable value of 
0.23. The distribution has the shape of a nearly symmetric, 
wide bump with a peak at Rw = 0.4. This is the simplest 
(single-parameter-based) way to put the stream model into 
consistence with data at a considerable probability level. 


5.3 Emission geometry plus selection effects 

An excellent agreement of the stream model with the data 
can easily be achieved when both the geometric parameters 
and selection effects are simultaneously taken into account. 
Fig. [8ji presents a case with streams located closer to the 
main meridian A <j>m = 0.4 x 90° = 36° and the blending 
unresolved below 0.2W O ut- Note the change of the Q-to-M 
profile fraction (cf. Fig. [TJ and the consistency of the total 
Rw distribution with data (KS probability of common ori¬ 
gin: 0.77). Fig. [SJr presents a similar result (A 4>m = 0.4x90°) 
for a fixed longitude resolution A(j> Tes = 1.5°. The probabil¬ 
ity of consistency with Fig. has increased up to 0.998. 
Note the peaks of the Q and M-type subdistributions in 
Fig- [Ur are narrow and misaligned which produces a double- 
peaked histogram. In Fig. [H)o the subparts are misaligned 
but wider, hence the shoulder at Rw = 0.3. 
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Figure 8. Simultaenous effects of component blending and az¬ 
imuthal confinement of streams. In both panels A (pm = 0.4 X 90°, 
ie. the streams are limited to two zones centered at the main 
meridian: cp m € (—36°, 36°) and (144°, 216°). a) A<p Ies = 
0.2Wout, b) A(p res = 1.5°. Note that such circum-meridional con¬ 
finement of streams increases the fraction of M type profiles (grey 
part of the histograms). 


6 CONCLUSIONS 

We find that the markedly different types of radio emission 
region (conal versus fan beam) imply dramatically different 
distributions of peak separation ratio for profiles with 4 and 
5 components. In the conal model, the ratio Rw is vastly 
dominated by values close to the size ratio R p of the two 
cones. This has been shown to be around 0.75 in several 
published analyses. In the case of the stream model, a much 
wider distribution with a peak at Rw ~ 0.2 is predicted. The 
observed distribution is a broad bump centered at Rw ~ 0.4. 

When these simplest, raw predictions are compared to 
the observed distribution, the stream-based model is many 
orders of magnitude more probable than the conal one. How¬ 
ever, this large difference between the predictions is mostly 
lost, when selection effects and geometric parameters are 
allowed to enter or change. The conal model can be made 
consistent with data by invoking our incapability to resolve 
components located closer than 0.22W / ou t■ Such resolution 
limit scales proportionally to the profile width, which could 
only be interpeted as an intrinsic effect of thick-walled emis¬ 
sion rings occupying a fixed fraction of polar tube. Such in¬ 
terpretation requires that a four times larger number of pul¬ 
sars (than the summed number of known Q and M profiles) 
should be hiding the merged pairs of components from our 
view. Such unresolved Q and M profiles should be observed 
as double and triple ones. 
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The stream model (fan-beam model) can be made con¬ 
sistent with data when the streams are positioned closer to 
the main meridian. A good agreement has been achieved 
for a zone consisting of two parts (poleward and equator- 
ward) centered at the main meridian, each of full width 
A0m = 72°. When the selection caused by the component 
blending is also taken into account, a perfect agreement of 
the stream model with data can be reached with no diffi¬ 
culty. 

The one-dimensional distribution of the peak separa¬ 
tion ratio is a new and interesting probing tool, which is 
free from (or mostly insensitive to) several unknown param¬ 
eters, such as emission altitude of emission region or dipole 
tilt. Nevertheless, the method depends on the observational 
capability to resolve blended components, which makes it 
less definite than the raw predictions of Fig. [I] Still, our 
results corroborate the success of the fan beam model in re¬ 
producing the properties of pulsar profiles. On the contrary, 
the conal model is allowed to persist only under specially 
tuned circumstances. 

The poor performance of the nested cone model under 
the Rw test, deserves a discussion in view of its previous 
success in reproducing the pulse-width distributions. The 
conal methods, as well as the recent method of Rookyard et 
al. (2015), assume that the lateral width of the radio beam 
decreases with the increasing impact angle /3, which is not 
the case for the fan beam model. The method of Rookyard 
et al. allows for a beam which is to a large degree arbi¬ 
trary (eg. patchy), however, since they associate the leading- 
or trailing-side outskirts of a pulse profile with the circular 
crossection of the polar tube, their method shares important 
qualitative properties with the conal model. All those meth¬ 
ods imply a nonisotropic distribution of the dipole tilt, with 
moderately small values of a preferred (Tauris & Manchester 
1998; Rookyard et al. 2015; and the references in Tab. [T]). 
The anisotropic a distribution is a valid possibility, since the 
radio emissivity may depend on a, just as expected from the 
accelerating electric field (Arons 1983; Harding & Muslimov 
1998). 

Our method assumes the isotropic a distribution, how¬ 
ever, this assumption may not be biasing our results, because 
the Rw method is essentially independent of a. It is there¬ 
fore possible that the small-a preference, as found by the 
other methods, results from the incorrect (circular) shape of 
the beam’s outer boundary. Qualitatively, this explanation 
would work in the right direction, for the following reason. 
Profiles produced by fan-shaped beams usually widen with 
the increasing impact angle /3. This is the case when a single 
outflowing stream extends laterally over an interval of the 
magnetic azimuth, or when there is more than one stream. 
For a central passage of the sightline (/3 = 0) the width of a 
profile formally vanishes (W = 0). Therefore, to reproduce a 

4 Detailed shape of the Rw distribution depends on a when the 
components’ resolvability is limited by the fixed A0 res (Sec. 11.11 . 
The bump at a moderate Rw is then present (Fig. [5j)), which 
corresponds to the cases with large a and f. Therefore, when the 
a distribution is assumed to concentrate at small values, where 
the selection effects cease to operate, the bump in the Rw dis¬ 
tribution becomes less pronounced. The small-a preference then 
makes the histogram more discordant with the observations than 
in the case of the isotropic a distribution. 


given observed pulse width W, a fan beam must be traversed 
at an appreciably large /3. This, through the observationally- 
fixed slope of a polarisation angle curve (S = sin a/ sin/3), 
implies larger a than for the conal beam. The conal beam 
model, on the other hand, may tend to underestimate the 
values of a and /3, in its effort to fit the observed widths of 
profiles (or core components) through the sightline travers¬ 
ing too close to the dipole axis. Unless we allow for the 
strong selection effects, the Rw method is only sensitive to 
parameters that determine the beam shape (not the scale). 
Therefore, we suggest the problems of the traditional conal 
model are inherent to the model itself, rather than to our 
method. The problems are probably caused by an incorrect, 
or at least not universal beam shape. 

A population of the nested-cone beams with a fixed R p , 
stands out as a narrow spike at Rw = Rp in the Rw distri¬ 
bution. The latter is observed to have the boxy shape with 
the Rw in the range between ~0.3 and 0.65. Therefore, the 
nested-cone beam can be made consistent with Fig. [Tjc, if a 
similar range of R p is assumed to exist in the real population 
of beams (with different values of R p being comparably nu¬ 
merous). However, this would imply an average R p of ~0.5, 
inconsistent with the numbers cited in Tab. [l] 

The present study has been inspired by the very sym¬ 
metric four-component profile of J0631+1036 (Zepka et 
al. 1996; Teixeira et al. 2016; Weltevrede et al. 2010). The 
inner components of its profile are located very close to each 
other, implying Rw = 0.33. As can be seen in Fig. 0J:, we 
have identified three more objects with Rw in the range be¬ 
tween 0.3 and 0.35. This makes up for 13% all objects in the 
observed Rw histogram. According to the conal Rw distri¬ 
bution (Fig.[4ji) only 1.6% of all Q and M objects should fall 
within 0.3 < Rw < 0.35. As noted in Teixeira et al. (2016), 
this suggests that the stream-based geometry (a system of 
fan beams) is responsible for the unusually symmetric profile 
of J0631+1036. This idea is strongly supported by the pres¬ 
ence of deep minimum at the center of the profile, (between 
the inner components), where the radio flux drops nearly to 
zero. This is difficult to understand within the conal model, 
because for Rw = 0.33 the nested cones with R p = 0.75 
imply the impact angle /3 = 0.95pi n , ie. at the center of the 
profile the sightline stays very close to the peak emissivity of 
the inner cone. Therefore, the nearly vanishing central flux 
seems to be incompatible with the standard conal geometry. 

Overall, the results discussed above provide more sup¬ 
port for the fan beam geometry, and raise more problems 
for the conal model. In view of the other arguments for 
the stream model (DRD10; Dyks & Rudak 2012; Wang et 
al. 2014; Chen & Wang 2014; Dyks & Rudak 2015) it may 
be worth to shift the current paradigm of geometric studies 
of pulsars from the conal structure towards the azimuthal 
arrangement of fan beams. 
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APPENDIX A: THE CONAL 

Rvc-DISTRIBUTION IN THE FLAT GEOMETRY 


When the beam is narrow (p ou t 1 rad) and the dipole tilt 
is large (a pout) the Rw distribution is well approximated 
by the flat geometry of Fig. [T] In that case the components 
are detected at the pulse longitudes: 

4>l n = Pin ~ P 2 , 4>lnt = Pout ~ p 2 , (Al) 

which determine the widths Win = 2</>i n and Wout = 2<^> ou t - 
The values of Rw = Win/W 0 nt = <(>in/</>out and the impact 
angle P are then given by: 


,2 _ Pin — P 2 o2 _ P?n(l “ R-W^p 2 ) 

w ~pL^P 2 ’ p “ 1 -R 2 w 


(A2) 


The number of observers that record Rw within some inter¬ 
val of A Rw is proportional to the interval of impact angle 
A/3 which corresponds to that A Rw- Therefore, 


n(Rw) 


lim n(ARw) oc 

AR\y -AO 


dp 

dRw 7 


(A3) 


which gives eq. ©. 


APPENDIX B: DEPENDENCE OF THE CONAL 
REDISTRIBUTION ON THE DIPOLE TILT a 


To simplify the calculation, we consider a central-cut case 
with £ = a (ie. /3 = 0). Pulse longitudes of the inner and 
outer pair of components are then given by: 


COS (j)i 


cos pi — cos 2 a 
sin 2 a 


, cos pi - 1 

J- I ■} ? 

sin a 


(Bl) 


where i = ‘in’ or ‘out’, When eq. 03 is Taylor-expandend 
to the 2nd power of <fii and pi, one obtains 4>i ~ Pit sin a, so 
that Rw = Win/Wout = pin/pout = R p does not depend on 
a. This explains the stability of the Rw distribution visible 
in Fig. 0i. 

To recognize the weak a-dependence of Rw, the cos (pi 
on the left-hand side of eq. 03 needs to be expanded to 
the order of rj>f. This leads to the quadratic equation for <j>f\ 


<t>i J.2 , 2(1 - COS Pi) _ n 

1 o (Pi ' -2 ^ 

12 sin a 


(B2) 


which has the solution: 
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6 (l - V 1 ~ x i) 


1/2 


with xl = 2(1 — cos pi )/(3 sin 2 q) 
ingly: 


(B3) 

Pi/( 3sin 2 a). Accord- 


1 - \/l -x? n _ 1 - y/1 - x 2 out R 2 P 

1 - a/1 ~ x Lt 1 - a/ 1 ~ X °ut 


(B4) 


Expanding the square roots around x^ ut -C 1 (up to the 
terms oc xZnt) one arrives at 


Rt 


r: 1 - 


1 -Ri 


1 + 12 sin 2 a/p 2 


(B5) 


Neglecting the unity in the denominator, taking a square 
root and Taylor-expanding it up to p 2 ut / sin 2 a gives the 
approximate eq. ©• This result shows that the value of Rw 
depends on a only for a nearly-aligned geometry (a < p out ). 
Therefore, the Rw distribution is insensitive to the assumed 
distribution of a, unless the latter is extremely non-isotropic, 
eg. with most objects having a of a few degrees. 


APPENDIX C: PULSE LONGITUDES FOR 
EMISSION FROM A FIXED MAGNETIC 
AZIMUTH 

To compute Win, W ou t, and R-w for the stream model (fan 
beam model) one needs a prescription for how to calculate 
pulse longitudes (p that correspond to the points where the 
selected (radio-emitting) magnetic azimuths (pm are sampled 
by the line of sight. A simple way to do this is to use eq. (19) 
of DRD10: 

cos ( = cos(-7r — (pm) sin a sin 9 m + cos a cos 9 m (Cl) 

to calculate the polar angles dm between the dipole axis 
and the emission direction from the sampled points (ie. the 
points at which the tangentially-emitting streams are de¬ 
tectable by the line of sight). Then the pulse longitudes (p 
for each component (ie. for each crossing point with some 
magnetic azimuth c p m ) can be found in the usual way: 

cos (p = (cos 9m — cosacosC)(sinasin<e) -1 . (C2) 

This would have been the full procedure, had it not been for 
a few, following technical complications. 

First, the cosine theorem of eq. m leads to the follow¬ 
ing quadratic equation for cos 9 m - A cos 2 8 m +B cos9 m +C = 
0, where: 

A = cos 2 a + cos 2 (pm sin 2 a (C3) 

B = — 2 cos a cos £ (C4) 

C = cos 2 ( — cos 2 (pm sin 2 a. (C5) 

For the positive discriminant A, ie. for sin 2 £ > (1 — 
cos 2 (pm) sin 2 a, there are two real solutions for 6 m which 
are measured from the same magnetic pole. Since the re¬ 
sult depends on cos 2 (pm (eqs. IC3I IC5I) there is a fourfold 
degeneracy associated with (p m , ie. the cases with ±</> m or 
7r i (pm cannot be discerned. Eqs. m and (IC2l) are thus 
blind to the sign of the magnetic azimuth (ie. to the location 
of the stream on the leading or trailing side of the magneto¬ 
sphere) so the sign of (p must be manually set equal to the 
sign of the corresponding (p m (with the latter understood 


in the range of ±7r). Moreover, the result is insensitive to 
the continuation of a given magnetic meridian to the other 
side of the magnetic pole (the solution for (p m is the same 
as for (pm + 7r, ie. it does not depend on the sign of 9m)- 
For example, consider a stream that starts on the poleward 
side of the dipole axis and extends away and upwards, into 
the rotationally-circumpolar regions of the magnetosphere, 
at some fixed magnetic azimuth (pm- When the line of sight 
has sin £ > sin a, it is passing on the equatorward side of the 
dipole axis, so the poleward stream should be missed near 
the dipole axis, with no achievable solution for <p. However, 
the l (pm plus pi’ denegeracy will result in a real solution for 
the extention of the azimuth < p m to the other (equatorward) 
side of the magnetic pole ((pm + 7t). A simple way to reject 
these false solutions is to insert all the calculated 9 m into 
eq. ED to check if the implied values of (/test are consistent 
with the original (. 

For large dipole inclinations, the larger of the two solu¬ 
tions for cos 9 m (hence related to a small 9 m ) corresponds 
to the stream crossing at a ‘near’ magnetic pole. The other 
solution has a smaller value of cos 9 m , with 9 m always mea¬ 
sured from the same (near) magnetic pole. In this latter 
case one may have cos 9 m ~ — 1 and 9 m ~ 180° which corre¬ 
sponds to the passage through the same magnetic azimuth 
close to the other (far) magnetic pole. Since the radio emis¬ 
sion is assumed to be latitudinally-limited to a single^] cir¬ 
cumpolar range of < 9 m < 0/ILx •C 1 rad, the far 
solution for cos 9 m (the one for which 9 m ~ 7r) is usually re¬ 
jected. However, in the case of the nearly aligned geometry 
a ~ £ £ both the solutions for 9 m (the smaller and 

the larger one) can survive at a single magnetic pole, ie. our 
line of sight can sample each (p m twice, while staying within 
the radio-emitting zone of 9 m - 

Although we have had few nearly aligned (very wide) 
profiles in our observed sample (eg. B1831—04, see Tabled, 
the complications resulting from the nearly aligned geome¬ 
try (small a) have been carefully treated. With five streams 
extending both into the poleward and equatorward part of 
the magnetosphere, it is possible to record components with 
more than five streams in the nearly aligned geometry. Ten- 
component profiles are possible when the sightline crosses 
each stream twice, as well as in the case when a -C 9ma* an d 
all ten streams are crossed once per period. In such cases, the 
definition of the outer and inner pair components becomes 
somewhat arbitrary, so we assume that the off-pulse region 
encompasses the azimuth (p m = <p = 180°. The correct or¬ 
dering of components is then achieved when the leadingmost 
component simultaneously has the smallest pulse longitude 
and magnetic azimuth (when both are defined in the range 
(—180°, 180°)). In the case of a « 0, and a < ( & 0/Tax, 
the component with the smallest value of (p m is unique. In 
general, however, e.g. for a = 0Zd x /2 and ( < 0/(Lx/2, there 
may exist pairs of components with the same value of (p m , 
including two with the minimum (p m but different (p. 

Accordingly, our simulation of Rw, Win and Wout for 
the stream model has followed these steps: 1) Ten mag¬ 
netic azimuths have been selected (five in the equatorward 


5 The inclusion of the far magnetic pole would just renormalise all 
the distributions by a factor of two. Since we ignore the question 
of interpulses, the second pole is neglected. 
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range of <j>m £ (—90°, 90°), the other five within (90°, 270°). 
2) The values of a and £ have been isotropically selected; 
the cases with |/3| > have been ignored. 3) The two 
groups of solutions for 8 m from eq. CD have been calcu¬ 
lated; those which implied inconsistent £ have been rejected. 
Those which fell outside the [9™^, 8^ x ) interval were also 
ignored. 4) The pulse longitudes have been calculated from 
(IC2l) with the sign of <f> set the same as for the corresponding 
sin< pm. Note that taking the sign of sine instead of that of 
(f>m redefines (transforms) the interval of </> m from the ini¬ 
tial (—90°, 270°) (useful for poleward/equatorward stream 
selection) into (—180°, 180°) (useful for the component or¬ 
dering with the off-pulse at </> m = <j) = 7r). 5) Both groups 
of solutions have been merged into a single one-dimensional 
array, then sorted and indexed increasingly. This produces 
a one-dimensional vector of solutions with the offpulse re¬ 
gion encompassing <j> = (f>m = —180°. The leftmost (lead- 
ingmost) solution has the smallest azimuths (both the mag¬ 
netic and the rotational one) when both are defined in the 
range of ±7r. 6) Only the cases with 4 or 5 solutions (cor¬ 
responding to Q and M profiles) have been considered. The 
widths Wout and Win have been calculated as the phase dif¬ 
ferences between appropriate components: W ou t = 4> n — <l> i, 
Win = 4>n -i — <)>2 where n = 4 or 5. The central solution 
(core) was only used to calculate separations between ad¬ 
jacent components A^ a dj, needed to consider the selection 
effects that result from the component blending. Otherwise 
the ‘core’ has been ignored. 
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